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Isoscalar meson spectroscopy from lattice QCD 
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We extract to high statistical precision an excited spectrum of single-particle isoscalar mesons 
using lattice QCD, including states of high spin and, for the first time, light exotic J^'~^ isoscalars. 
The use of a novel quark field construction has enabled us to overcome the long-standing challenge 
of efficiently including quark-annihilation contributions. Hidden-flavor mixing angles are extracted 
and while most states are found to be close to ideally flavor mixed, there are examples of large 
mixing in the pseudoscalar and axial sectors in line with experiment. The exotic J^*^ isoscalar 
states appear at a mass scale comparable to the exotic isovector states. 



Introduction: Mesons which have all flavor quantum 
numbers equal to zero, known as isoscalars, form a par- 
ticularly interesting sector of QCD, the gauge-field the- 
ory of quarks and gluons. They can be constructed from 
quark-antiquark pairs of any flavor or even from states 
of pure glue, and in general the QCD eigenstates will be 
superpositions of these possible basis states which mix 
through quark annihilation. By studying the spectrum 
and hidden-flavor content of isoscalar mesons we can in- 
fer a phenomenology of annihilation dynamics in QCD. 

Empirically it is found that many low-lying isoscalar 
mesons come in identifiable pairs with a strong prefer- 
ence for ideal flavor mixing; prominent examples include 
u)^<f> where the admixture of \ss) into the dominantly 
H- \dd)) w is estimated to be below 1%. Sim- 
ilarly the tensor mesons, /2(1270), /2(1525), where the 
lighter state decays almost entirely into tttt and the heav- 
ier into KK suggesting mixing of at most 1% There 
are notable exceptions to the dominance of ideal flavor 
mixing in the pseudoscalar and axial sectors. The rj and 
r]' mesons appear to be close to the octet and singlet 
eigenstates of SU{3)p which have lowest quark content 
^ {\uu) + \dd) - 2\ss)) and ^ {\uu) + \dd} + \ss)) re- 
spectively and are thus far from ideally flavor mixed. 
The scalar sector poses the greatest challenge to phe- 
nomenology - the large number of low-lying resonances 
has suggested interpretations featuring light and strange 
QQ^ QQ'l'l states, glueballs and meson-meson molecules ■ 

Our aim here is to extract an isoscalar meson mass 
spectrum and study the hidden-flavor composition of 
the states using lattice QCD computations. This is a 
challenging undertaking for several reasons. It requires 
lattice gauge configurations with dynamical light and 
strange quarks in order that the flavor mixing appear in 
a unitary manner, while evaluation of the disconnected 
correlator contributions that distinguish isoscalars from 
isovectors is computationally expensive and the signals 
obtained typically diminish into noise at small Euclidean 
times. These problems have limited calculations to a few 



with typically low statistical precision Glue- 
ball studies, which produce exceptionally clean spectra 
in the quark-less Yang-Mills case 16| , become challenging 
in QCD through strong coupling to the quark sector Q. 

In this letter we present the isoscalar meson spectrum 
for a single choice of light and strange quark masses, lat- 
tice spacing and volume across multiple J^^ , excluding 
the 0"*""*" case which is of sufficient interest to justify a 
separate publication of its own. This extends the work 
reported in 0, for isovector and kaonic mesons, taking 
advantage of many of the techniques developed therein. 

Lattice technology: We compute the spectrum of 
isoscalar mesons using a basis of operators of the form 
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d, and s are the up, down and strange quark fields, and 
the rf are operators acting in space, color, and Dirac 
spin space [10| on a time slice, t. We combine these 
operators to construct two-point correlators of the form 
C%l^{t',t) ^ {0\O\{t')O'j^{t)\0) which, after integration 
of the quark fields, can be composed from connected com- 
ponents, diagonal in quark fiavor, 

and disconnected components that can mix flavor, 

viUt'^t) = Tr[<f>^(t')v(t',t')] Tr[$^(t)r,(t,0]. 

The quark flclds in O arc acted upon by a "distilla- 
tion" smearing operator that emphasizes the low mo- 
mentum quark and gluon modes that dominate low mass 
hadrons 0, This smearing can be factorized al- 

lowing the "perambulators", r^, and the zero momen- 
tum projected operators, $, to be constructed as matri- 
ces in distillation space: Tq{t',t) = Vt]M^\t',t)Vt and 

<I>^(t) = V^TfVt, where Vt are the eigenvectors of the 
gauge-covariant three-dimensional Laplacian operator on 
a time-slice t, and Mg is the lattice representation of the 
Dirac operator for quark flavor q. The resulting traces are 
then over the set of eigenvectors which is much smaller 



than the full lattice space - this allows for an efficient 
computation of the disconnected terms, T), that is not 
possible in the traditional "point-to-all" method. 

We consider the limit where the u and d quark masses 
are set equal, and are replaced by a "light" quark, £. 
For a particular combination of creation and annihilation 
operators {A,B), the flavor subspace in the correlation 
matrix is then 
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The signals for V and C are typically short in physical 
units. To help resolve this, wc use anisotropic lattices 
with three dynamical flavors of Clover fermions [111 [12 1 
of size 16"^ X 128, spatial lattice spacing Us ^ 0.12 fm, 
and temporal spacing, a^^ ^ 5.6 GcV. This fine tempo- 
ral lattice spacing has proven useful in determi ning the 
spectrum of isovector mesons [1,0] and baryons flij . 

The disconnected terms V require perambulators 
Tq{t,t) with sources on every time-slice. With Tq{t',to) 
from all tinieslice sources, to, in hand, we can maxi- 
mize the signal above noise by averaging: '^{t) = 

where the disconnected term is averaged over all = 128 
timeslices, while the connected term needs far less averag- 
ing to achieve a comparable level of statistical fluctuation 
- in this study we used A^t = 32. Since the eigenvectors 
Vt sample the full three-dimensional space and the in- 
versions are from all time-slices, the propagation from a 
source is from every point of the four-dimensional lat- 
tice, resulting in a signiflcant decrease in the noise of the 
correlator compared to "point-to-all" methods. In all, 
479 configurations (separated by 20 trajectories) were 
used with 64 distillation vectors. Graphics Processing 
Units 3, l3l allow us to evaluate the 31-million Dirac 
propagator inversions required An example of the 

correlator quality is shown in Figure [T] for the case of the 
distilled 75 operator. The disconnected contributions are 
resolved for time-separations beyond 1 fm. 

Extraction of the spectrum: Our approach is to form 
a matrix of correlators with elements C'^g{t) that is 
then diagonalized in the (operator (g) flavor) space to 
obtain a spectrum which best describes the correlators 
in the variational sense @. This yields both estimates 
of the mass spectrum, {m„}, and the overlap of states 
onto the operators used to interpolate from the vacuum, 
Z^^^ = (n|O^'^(0)|0). The overlaps can then be studied 
to explore the hidden-flavor composition of the states. 
We use an operator set founded upon the derivative- 
based constructions for F described in Q, doubled in 
size by forming operators from both light-quark {q = £) 
and strange-quark (g = s) fields. 

The variational analysis is performed independently in 
each cubic-group irreducible representation, A^*^, result- 
ing in spectra that, as in [§], conform to the mass de- 
generacy and overlap patterns expected of single-particle 
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FIG. 1: Time-source averaged connected and disconnected 
correlators with Fa = Fs = 7^. 
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FIG. 2: Effective mass of lowest two variational principal cor- 
relators in the pseudoscalar channel, m]^g = — In |^-3jjy-^j , 
with St — 3. The non-monotonic behavior of the statistical 
error is an artifact of our implementation of the variational 
analysis method. The horizontal bands indicate the mass re- 
sult from fits to the principal correlators. 



states of good J^'-^ subduced into the relevant A^'-^. We 
find that while isoscalar signals arc typically noisier than 
the corresponding isovector correlators, they arc not ex- 
cessively so and a detailed spectrum comparable to that 
presented in Q for isovector states follows. As an exam- 
ple of the quality of spectrum signals observed, in Figure 
[2|we show the effective mass of the principal correlators 
for the lightest two states with ^ quantum numbers 
which wc identify with the ?/ and ry' mesons. 

Within a given J^^ sector of the extracted spectrum, 
we find that states occur in pairs with overlaps that 
suggest they are admixtures of the light-quark, strange- 
quark basis states. Wc parameterize this admixture by 
introducing a mixing angle, assuming the pair of states 
(a, b) are orthogonal combinations of just two light- 
strange basis states: 



where \£) = {^\uu) + \dd)) , \s) = \ss) and where 



a) = cos a 


— sin a 


b) = sin a 


t) + cos a 
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FIG. 3: Effective mixing angle, aA{t), as a function of Eu- 
clidean time for the lightest two pseudoscalar states. Colors 
indicate seven different operator constructions, Fa, and the 
grey band a constant fit. 



a is the lighter of the two states 
is determined 12: 



The mixing angle 

A 



tan 



for each operator, T , from = 
where the overlaps Z are deter- 



mined on each timeslice by the variational diagonaliza- 
tion and where we're assuming that e.g. (^lO^'^lO) = to 
a good approximation. In general we find a good agree- 
ment between determinations from the different opera- 
tors forming our variational basis - see Figure [3] for the 
example of the 77, rj' states. 

The quark-bilinear operator basis thus far described 
is imperfect for the study of isoscalar mesons since it 
lacks constructions likely to have a large overlap onto 
glueball basis states. Such operators can be constructed 
from gauge-invariant link paths on the lattice, and it 
is known that in the quark-less Yang-Mills S'C/(3)-color 
theory there is a spectrum of bound states, the lightest 
of which fall within the mass range considered in our 
study [1] . We included operators of this type in the basis 
but it was found that correlators featuring them had sig- 
nals which degraded into noise at rather short Euclidean 
times before any clear state signals could be extracted. 
Any resolution of such states will require operators that 
plateau at shorter Euclidean times as well as more gauge 
configurations. The presence of torelons, states of gluonic 
flux wrapped around the periodic spatial volume, while 
anticipated in the studied mass region, arc not observed 
in the spectrum. 

In Figure S] we present the main result of this study, 
the spectrum of isoscalar mesons with m^r ~ 396 MeV on 
a ~ (2fm)^ lattice identified by J^^ quantum numbers. 
The degree of light-strange mixing is indicated by the 
black-green coloring and by the displayed mixing angle 
a. Superimposed on the plot is the spectrum of isovector 
mesons on the same lattice, taken from reference and 
for comparison the spectrum of glucballs in the quark-lcss 



Yang-Mills theory taken from [6[. 

Phenomenology: The m,r ~ 400 MeV spectrum re- 
sults presented in Figure 2] suggest a phenomenology of 
isoscalar mesons that is in quite good agreement with ex- 
periment. We reproduce the pattern of rj, rj' observed in 
experiment and suggested by the axial anomaly of QCD, 
where the lighter 77 meson is dominantly an SU{3) p octet 
and the significantly heavier rj' is dominantly SU{?i)f 
singlet. In fact our estimate of the light-strange mixing 
angle, 42(1)°, deviates somewhat from the exact SU{3)f 
expectation of 54.7°, in line with phenomenological de- 
terminations (e.g. |16l4l8| and references therein). 

The level of statistical precision in the disconnected 
contribution to the vector channel is such that we can 
determine the uj — p mass splitting to be 21(5) McV. The 
very small 1.7(2)° uj/cp mixing angle is comparable to the 
3.2(1)° inferred from analysis of radiative decay rates of 
these states to pseudoscalar mesons [l8| . 

We observe that the bulk of states are close to ide- 
ally flavor mixed, and we find that the majority of 
light isoscalar states are heavier than their corresponding 
isovector cousins even in the cases of negligible admix- 
ture of heavier ss components. This is of course possible 
through the disconnected contribution to the correla- 
tor and hence the effect of "annihilation dynamics" . An 
exception to this observation is the 1'' channel, where 
the lightest hi, with negligible ss component, is lighter 
than the isovector 61, in line with experiment, although 
given the large hadronic decay width of the hi this may 
not be significant. Apart from the previously discussed 
pseudoscalars, the only other case of large mixing an- 
gles in a conventional J^*^ channel are the axial (l"*"*") 
states. For the lightest two states we find a mixing an- 
gle of 31(2)° which is in reasonable agreement with the 
phenomenologically determined 21(5)° between the 
/i(1285) and the /i(1420). 

The rightmost section of Figure |4] shows our determi- 
nation of the spectrum of exotic J^'~^ isoscalar mesons. 
This is the first computation of such a spectrum and indi- 
cates that isoscalar exotics appear at a comparable mass 
scale to the exotic isovectors, and as in @, the pattern 
of operator overlaps suggests that these states are hybrid 
mesons. While the 0"' and 2"' states are rather close 
to ideally flavor mixed, there is considerable mixing in 
the 1 ^ channel. 

Outlook: In summary, we find that lattice QCD can 
reproduce the pattern of single-particle light isoscalar 
states observed in nature, and the determined flavor mix- 
ings are in reasonable agreement with phenomenology. 
Beyond their role in the isoscalar spectrum, annihilation 
dynamics feature prominently in hadron decays. Com- 
bining the methods developed in this letter with finite- 
volume techniques for the extraction of phase shifts [l3| 
and stochastic methods [i^l, future work will concen- 
trate on the determination of hadronic resonances within 
QCD. 
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FIG. 4: Isoscalar meson spectrum labeled by J . The box height indicates the one sigma statistical uncertainty above and 
below the central value. The light-strange content of each state (cos^ a, sin^ a) is given by the fraction of (black, green) and 
the mixing angle for identified pairs is also shown. Horizontal square braces with ellipses indicate that additional states were 
extracted in this J^'^ but were not robust. Grey boxes indicate the positions of isovector meson states extracted on the same 
lattice (taken from [^). The mass scale is set using the procedure outlined in [3, Pink boxes indicate the position of 

glueballs in the quark- less Yang-Mills theory [^. 
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